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This study combines surface-sensitive photoemission experiments with density functional theory 
(DFT) to give a microscopic description of H adsorption-induced modifications of the ZnO(lOlO) 
surface electronic structure. We find a complex adsorption behavior caused by a strong coverage 
dependence of the H adsorption energies: Initially, 0-H bond formation is energetically favorable 
and H acting as an electron donor leads to the formation of a charge accumulation layer and to 
surface metallization. The increase of the number of 0-H bonds leads to a reversal in adsorption 
energies such that Zn-H bonds become favored at sites close to existing 0-H bonds, which results 
in a gradual extenuation of the metallization. The corresponding surface potential changes are 
localized within a few nanometers both laterally and normal to the surface. This localized character 
is experimentally corroborated by using sub-surface bound excitons at the ZnO(lOlO) surface as a 
local probe. The pronounced and comparably localized effect of small amounts of hydrogen at this 
surface strongly suggests metallic character of ZnO surfaces under technologically relevant conditions 
and may, thus, be of high importance for energy level alignment at ZnO-based junctions in general. 

PACS numbers: 73.20.At, 73.20.Hb, 79.60.Dp, 73.61.Ga 


I. INTRODUCTION 

Zinc oxide (ZnO) is a wide band gap {Eg = 3.4 eV), 
intrinsically n-type, transparent conductive oxide that 
has received much attention due to its potential use 
in novel opto-electronic devices such as organic light- 
emitting diodes and photovoltaics.^®^ Moreover, it is al¬ 
ready widely used in catalysis and chemical sensingIn 
all these cases, device functionality is governed mainly by 
the properties of the ZnO surface or its interface with, 
e.g., functional organic molecules. The opto-electronic 
properties of such interfaces depend critically on the 
alignment of energy level^ and the occurrence of collec¬ 
tive surface phenomena such as ZnO surface excitons 
ZnO also has a tendency to grow in self-organized 
nanoscale structures (rod s, wir es, ribbons, etc.) with 
large surface-to-bulk ratios.^^SE^In these nano-structures, 
the mixed-terminated (1010) surface dominates, because 
it is energetically the most favorable.tl^ The optical, elec¬ 
tronic and catalytic properties of ZnO are highly sensi¬ 
tive to surface modifications such as impurities, defects 
or adsorbates,^ and yet a fundamental understanding of 
the key phenomena at the ZnO surface remains elusive.^ 
This gap in our understanding is also a problem for sur¬ 
face functionalization by attachment of optically active 
molecules, where knowledge about the interfacial elec¬ 
tronic structure is often restricted to a macroscopic, av¬ 
eraged view. In this Article we address this knowledge 
gap, by focusing on hydrogen (H) adsorption and by com¬ 
bining surface-sensitive photoemission experiments with 
ab initio density functional theory (DFT) to give a con¬ 
cise microscopic view on H adsorption-induced effects at 
the ZnO(lOlO) surface. 

In recent years, the significant impact of H doping 
on the electronic structure of ZnO has been recognized. 


which is highly relevant since hydrogen is a ubiquitous 
and prominent contaminant in both laboratory and in¬ 
dustrial settings. Using density functional theory (DFT) 
it was found that atomic H acts as a source of the (unin¬ 
tentional) n-type condu ctivity by forming shallow donor 
states in bulk ZnO.This prediction was corroborated 
by several studies and extended from interstitial to ad¬ 
sorbed hydrogen at the ZnO surface.t^^lli^ 

The donor character of adsorbed H on ZnO has been 
known for decades. The experimental observation of in¬ 
creased conductance of ZnO films due to H adsorptiorP^ 
has been exploited by using ZnO as h igh-se nsitivity sen¬ 
sor for hydrogen and carbon hydratesIn the estab¬ 
lished picture, the electron donor character of H induces 
downward surface band bending that leads to a crossing 
of the conduction band minimum (CBM) and the Fermi 
level As a consequence, a strongly confined surface 
charge accumulation layer (CAL) forms,l^as sketched in 
Fig. a). At the same time, this adsorption-induced in¬ 
terfacial charge transfer leads to a reduction of the net 
surface dipole, which results in a reduction of the work 
function 4>. Moreover, analogous H-induced CAL forma¬ 
tion was found in other transparent conductive oxides, 
SrTiO3(001), SnO2(110), and MgOPEH in contrast to 
conventional semiconductors, e.g.. Si or GaAs, the sur¬ 
face band bending on n-type ZnO extends only a few 
10 A into the bulk, making it a exclusively surface-related 
phenomenon.!^ First direct spectroscopic evidence of H- 
induced CAL formation at the ZnO(lQlQ) surface was 
shown by Ozawa and Mas (j^^ * ^^ who observed occupied 
states at the F-point using angle-resolved photoelectron 
spectroscopy (ARPES). The authors proposed that the 
surface metallization is a consequence of the partial oc¬ 
cupation of Zn45 conduction band states. This results 
in a single metallic band confined in the potential well 
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between the vacuum interface and the CBM, in which 
a laterally delocalized two-dimensional electron gas with 
a maximum charge density of 10^^ cm“^ is formed [cf. 
Fig. a)]. It was shown by Woll and coworkers, that the 
formation of 0-H-bonds at the (1010) surface is respon¬ 
sible for this semiconductor-to-metal transition.^^ Based 
on their calculations they conclude further, that a com¬ 
plete saturation of both 0-H and Zn-H-bonds at the 
surface leads to the restoration of the semiconducting 
state. 



FIG. 1. (a) Schematic energy level diagram of ZnO(lOlO) 

after exposure to H with an average potential V(z) across the 
interface. The CBM, which is located only 200 meV above 
the Fermi level in n-type ZnO(lOlO), is bent below Ef leading 
to the formation of a few 10 A thick CAL. The net surface 
dipole is modified by H adsorption. Above-band-gap pho¬ 
toexcitation creates surface excitons that can be measured by 
photoemission with a time-delayed second pulse as indicated 
by the arrows, (b) Representative photoemission spectrum of 
H/ZnO(1010) showing the secondary electron cutoff at 
Fkin = 0 and the occupation of electronic states in the CAL 
below Ef- 

Previous ZnO(lOlO) studies have focused on surfaces 
with H adsorbed either exclusively on surface O sitesor 
have used comparably high H dosages that lead t o a sat- 
uration of all energetically possible surface sites,^^-^^^^!^ 
These studies do therefore not investigate the competi¬ 
tion of H adsorption at the two different sites, especially 
in the regime of low coverages. Although the different 


effects of 0-H vs. Zn-H bond formation have been rec¬ 
ognized previously, a detailed microscopic view of how 
single H atoms interact with the pristine or slightly H- 
covered ZnO(lOlO) surface is not established. The pic¬ 
ture of a laterally delocalized CAL can only provide a 
macroscopic, averaged description of the surface poten¬ 
tial changes [see Fig. [^a)]. 

In the present study, we focus on the low H cover¬ 
age regime to elucidate the impact of H adsorption on 
the local potential and electronic structure of the mixed- 
terminated ZnO(lOlO) surface. We measure the elec¬ 
tronic states of single crystal ZnO(lOlO) surfaces using 
laser-based photoelectron spectroscopy (PES). This in¬ 
herently surface sensitive method gives direct access to 
the occupied electronic states and work function of the 
ZnO surface. Complemetarily, we employ ab initio DFT 
calculations which facilitate an atomic scale description 
of the electronic structure. We demonstrate that differ¬ 
ent H pre-coverages, i.e., the amount of H already ad¬ 
sorbed on the surface, lead to drastic changes in the ad¬ 
sorption behavior of further H atoms. Furthermore, we 
develop a microscopic picture of the markedly different 
effects that 0-H and Zn-H bond formation has on the 
landscape of the surface potential. Experimentally, the 
localized character of the H-induced potential changes is 
confirmed by using the sub-surface bound exciton as a 
local probe for the potential. These findings extend the 
conventional view of the delocalized and uniform charac¬ 
ter of the H-induced CAL at the ZnO(1010)-surface to 
the low-coverage regime. 


II. METHODS 
A. Experimental Methods 

The experiments and sample preparation were per¬ 
formed in an ultrahigh vacuum (UHV) chamber with 
a base pressure below 1 x 10“^^ mbar. Hydrother- 
mally grown single crystal ZnO(1010) samples (MaTecK 
GmbH) were prepared by Ar^ sputtering (0.75 keV, 
10 min) and annealing cycles at 750-850 K for 30 min 
with comparably slow heating and cooling rates of 
20 K min“^, following established procedures.^ Hydro¬ 
gen was offered at a constant background pressure of 
6.7(2) X 10“^ mbar (0.5 L s“^), at a sample temperature 
of 100 K for all measurements presented in this paper.l^ 
Atomic H was generated using a glowing tungsten fila¬ 
ment at a distance of ~ 15 cm in the line of sight of the 
sample surface, partly cracking the H 2 molecules. This 
is analogous to the procedure described in Refs. [T7| and 

inpii 

Eemtosecond laser pulses were generated by a 200 kHz 
regeneratively amplified laser system. Eor static photoe¬ 
mission measurements photon energies of hn = 4.63 eV 
and 6.18 eV were created by frequency tripling and qua¬ 
drupling of the fundamental 1.55 eV laser output. Eor 
time-resolved pump-probe photoelectron spectroscopy 
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light pulses with 4.19 eV photon energy were generated 
by frequency doubling the output of an optical paramet¬ 
ric ampliher and were then used to excite electron-hole 
{e-h) pairs [see Fig. [^a)]. The temporal evolution of 
this nonequilibrium state is monitored by a time-delayed 
probe pulse with a photon energy of 4.63 eV (see Ref. |9] 
for further details). Photoelectrons were detected using 
a hemispherical photoelectron analyzer (PHOIBOS 100, 
Specs GmbH) with energy resolution of 30-50 meV. The 
photoelectron spectra were integrated over an angle of 
^ ±2° around the F-point. The binding energy of the 
photoelectrons was referenced to the Fermi level (^f) of 
the Ta sample holder which was in electrical contact with 
the sample surface. 

The quality of the pristine surface was routinely 
checked by measurements of the work function = 
4.50(5) eV and the valence band maximum (VBM) at 
E — Ep = —3.18(6) eV, ah agreeing well with literature 
values .1^ The work function was determined using PES 
by measuring the energetic position of the secondary elec¬ 
tron cutoff which is constituted by photoelectrons 
that barely overcome ^ = hu — (Fp — F^^) as depicted 
in Fig.j^b). 

While illuminating the freshly prepared surfaces, we 
always observe a work function reduction by few 10 meV 
along with a slight increase of the CAL intensity on a 
timescale of several 10 s. This is attributed to the for¬ 
mation of a small portion of surface defects by the UV 
illumination.!^ Experimental results shown here always 
refer to stabilized conditions. Additionally, there is a 
comparable shift of 4> on a timescale of hours due to the 
ubiquity of H 2 as residual gas even in UHV. This can be 
neglected in our experiments. We carefully checked for 
spectral shifts due to charging or photovoltage effects by 
varying the incident photon flux and found no shifts, nei¬ 
ther for clean, nor for the hydrogen-covered ZnO surfaces. 
Low energy electron diffraction measurements showed the 
expected rectangular (1 x l)-pattern. 


B. Computational Methods 

Ah calculations were performed with the Eritz-Haber- 
Insti tute ab initio molecular simulations (EHI-aims) 
code.l^^Hll Unless otherwise noted, the functional of 
Perdew, Burke, and ErnzerhoP^ (PBE) was employed. 
We account for van der Waals forces through the scheme 
of Tkatchenko and ScheffleJ^ using the parameteriz ation 
as described in Ref. EZl The ZnO (1010) surface was mod¬ 
eled with a 4 X 4 unit cell containing 16 ZnO surface 
dimers, with a depth of 32 single-layers. We found that 
such a slab thickness was necessary to capture the whole 
extent of downward band-bending in the most extreme 
case in which each surface oxygen is decorated by a hy¬ 
drogen atom. Since band-bending is entirely enclosed 
in our superceh, we did not employ the electro static 
schemes recently developed by some of A re¬ 

gion of 30 A vacuum was inserted between the ZnO slab 



FIG. 2. Dependency of the work function of different 
ZnO (1010) samples on hydrogen dosage (blue curve) and in¬ 
tegrated intensity of the CAL signature (red circles). 


and its periodic replica. Polarization through the vacuum 
was prevented by means of a dipole correction.!^ The self- 
consistent held cycle was converged to 10“^ eV for the to¬ 
tal energy, 10“^ eV for the electron density and 10“^ eV 
for the sum of eigenvalues. In FHI-aims, the basis is hier¬ 
archically ordered in tiers^ For Zn we employed the tier 
1 basis and for O the tier 2 basis together with tight de¬ 
fault settings. For integrations, a tightly converged Lebe¬ 
dev grid was used. Ah calculations were done assuming 
a constant ZnO doping concentration of ~ 10^^ e cm“^, 
which was modeled using the virtual crystal approxi¬ 
mation approach.SIl Using the multis cale virtual crystal 
approximation dopant approachwe tested carefully 
for ah systems reported here that no transfer of bulk 
charge carriers to the quantum-mechanicahy treated 32- 
layer thick slab occurs. For the limiting case of a full O- 
H monolayer, we carefully verihed that the amount and 
the spatial evolution of the band bending is not affected 
by the choice of the functional by comparing to hybrid 
density function al calc ulations using the Heyd-Scuseria- 
Ernzerhof (HSE).!^^^^ We adjust the amount of exact ex¬ 
change to 0.4, in accordance with previous worlP^EHand 
denote this functional HSE*. Moreover, we compared 
the relative energies of a slab with 100% 0-H coverage 
and 0% Zn-H coverage and a slab with 75% 0-H and 
25% Zn-H using PBE and HSE* for a 2 x 2-supercell. In 
PBE, exclusive adsorption on oxygen atoms is 0.16 eV 
more stable than in HSE*. Hence, when comparing ad¬ 
sorption on O and Zn, we expect PBE to err slightly in 
favor of 0-H adsorption. However, this does not affect 
the qualitative hndings discussed in this work. 


III. RESULTS AND DISCUSSION 

An exemplary photoemission spectrum of the 
ZnO(1010) surface after dosing 25 L of H 2 is shown in 
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Fig.^b). At this H 2 dosage, we observe a work function 
of 4.03(5) eV and a peak right below = 0 eV, which 
we attribute to the CAL. To quantify the H-induced 
changes to the surface potential, we took a series of 
spectra for H 2 dosages ranging from 0 L to 500 L 
and analyzed the shift in 4> and the intensity of the 
CAL signature, which is shown in Fig. The blue 
curve depicts the dependency of 4> on the H 2 dosage. 
The work function is continuously reduced by up to 
A4>max = —0.65(5) eV, as compared to the work 
function of the pristine surface after annealing. Clearly, 
the shift of 4> occurs mainly at hydrogen dosages below 
50 L, and for dosages approaching 100 L we observe 
a stabilization of the work function. A quantitatively 
identical dependency of 4> on H dosage is measured for 
different samples, as shown in Fig. Our measured 
value of A4>max is identical to the value for H-saturated 
surfaces reported in the liter at ure.^^^ 

To evaluate the effect of H adsorption on the surface 
electronic structure we analyzed the intensity of the CAL 
signature with respect to the hydrogen dosage as shown 
by the red squares in Fig. [2p^ While the PE spectrum of 
the freshly prepared surface shows nearly zero intensity 
below Ep, the CAL intensity shows a distinct increase for 
hydrogen dosages up to 150 L. It should be noted that 
the spectral shape of the peak and its energetic position 
do not significantly change in this low dosage regime.l^In 
addition, the increasing intensity of the CAL peak at Ep 
implies an increase of the electron density and, therewith, 
metallic character of the ZnO surface. 

Continued dosing of hydrogen above 150 L then leads 
to a significant intensity reduction of the CAL signature, 
which is in contrast to the stabilization of 4> for high 
dosages. Despite this difference in behavior for high cov¬ 
erages, the characteristic change of the two H-induced ef¬ 
fects in the dosage range between 50 and 150 L hints at a 
competition between two different processes: (i) the cre¬ 
ation of DOS below Ep, i.e., metallization and a simulta¬ 
neous reduction of the surface dipole for dosages ^ 150 L 
and (ii) the subsequent reduction of the CAL intensity 
accompanied by a stabilization of the surface dipole for 
higher dosages. These opposing effects may stem from 
the availability of two different adsorption sites: Zn and 
O surface atoms. 

To gain microscopic insight into the adsorption ener¬ 
getics, we performed DFT calculations. Strictly speak¬ 
ing, the correct way to determine the surface structures 
would be to generate all possible hydrogen configura¬ 
tions and to calculate their surface energy. The differ¬ 
ent configurations would then be populated according to 
Boltzmann-statistics at the experimental temperature of 
100 K. The large number of possible configurations [a 
(4 X 4)-supercell can represent on the order of 2 x 10^ 
different geometries], renders this approach intractable. 

The traditional approach is to deal with this is to ne¬ 
glect temperature effects, i.e., to assume T = 0 K. This 
would reduce the problem to finding the conformation 
with the lowest energy, which could be done by cluster ex- 


next H 
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on Zn on O 



FIG. 3. (4 X 4) ZnO supercell with 25% 0-H pre-coverage. 
In our calculations we compare the energy difference for the 
adsorption of an additional H atom on either an O atom or a 
Zn atom next to an H-pre-covered O atom, as indicated. 



FIG. 4. (a) Relative formation energy for the formation of 

one 0-H bond (AG^) and one Zn-H bond (AG^"") per 4x4 su¬ 
percell, depending on the amount of oxygen atoms pre-covered 
with hydrogen. Formation of 0-H bonds is clearly favored for 
pristine ZnO (1010). The crossing of the zero line marks the 
0-H-bond pre-coverage at which the next H favorably adsorbs 
at a surface Zn atom next to an 0-H bond. This coverage of 
18 % is marked with a red dotted line, (b) Galculated work 
function change resulting from H adsorption at a surface with 
the same 0-H pre-coverage as in (a). Diamonds and circles 
show A4> for H adsorption on O or a first Zn site, respectively. 


pansion or genetic algorithm methods. However, picking 
several random configurations we found that the energy 
differences are very small, on the order of 10 to 50 meV. 
Assuming T = 0 K also neglects configurational entropy, 
which is not constant for a coverage series. Instead, con¬ 
figurational entropy favors submonolayers with medium 
coverage over almost full or empty layers. Moreover, for 
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a given hydrogen decoration, it favors layers with mixed 
adsorption sites over layers with only one adsorption site. 
Given the small energy differences between the geome¬ 
tries, we chose to simplify the calculations by assuming 
that for a fixed number of 0-H and Zn-H bonds all con¬ 
figurations are essentially degenerate. This relieves us 
from the burden of finding the global minimum, and we 
can simply pick one arbitrary geometry as representa¬ 
tive. Moreover, we can then use an analytic expression 
for the configurational entropy S. According to Boltz¬ 
mann’s equation, it is given as 


S' = kB\n 


ff 
V no y V ^Zn y 


( 1 ) 


Here, no and nzn denote the number of surface O and 
Zn atoms, respectively, and non and nznH the number 
of 0-H and Zn-H bonds. 

To keep the computational effort reasonable, we refrain 
from performing all 256 calculations that would be re¬ 
quired to generate a full surface energy diagram as func¬ 
tion of 0-H and Zn-H coverage. Instead, we focus on 
the question whether at some point during the dosage 
of hydrogen atoms, the formation of Zn-H bonds could 
energetically compete with the formation of 0-H bonds. 
Previous work by Woll and co- worker^ suggested that 
surface 0-H bonds are more stable than Zn-H bonds, 
since the latter were not observed at high temperatures. 
Furthermore, the sticking coefficient of hydrogen on the 
pristine Zn-terminated ZnO(OOOl) surface was found to 
be extremely small at < 1 x 10“^ (Ref. [T7|) . which il¬ 
lustrates the low formation probability of isolated Zn-H 
bonds. Taking this into consideration, we calculated the 
Gibbs energy for the formation of adding another hydro¬ 
gen atom to a surface that is already pre-covered with a 
sub-monolayer of hydrogen exclusively adsorbed on oxy¬ 
gen. We considered pre-coverages in a wide range from 
0% to 94%, specifically one, two, three, four, six, and 
eight H atoms in the unit cell, which were always dis¬ 
tributed such that the distance between them was maxi¬ 
mized (cf. Figure [^. Then, an additional hydrogen was 
introduced either on a Zn atom next to an 0-H bond or 
on an additonal O atom. An example geometry for a pre¬ 
coverage of 25 % is shown in Fig. The Gibbs energy 
of formation, AG is calculated separately for adsorption 
on Zn or O site as 


A^Zn/o ^ - TS (2) 


where F^n+i is the energy of the 4x4 supercell with one 
additional hydrogen adsorbed on either Zn on O, the 
energy of the supercell without the additional hydrogen, 
£'(H) the energy of a hydrogen atom, T the temperature 
in experiment (100 K) and S the configurational entropy. 
Note that we have neglected the contribution from vibra¬ 
tional zero-point energies here.^^ 

We find that with increasing (pre) coverage of 0-H 
bonds, the formation of further 0-H bonds becomes in¬ 
creasingly energetically unfavorable with respect to the 


formation of Zn-H bonds. We attribute this to the am¬ 
photeric character of hydrogen which acts as an electron 
donor on oxygen but as an electron acceptor on Zn. Pre¬ 
covering the surface with 0-H bonds reduces the work- 
function of ZnO and thus also the ionization energy and 
electron affinity. Gonsequently, with increasing 0-H pre¬ 
coverage, charge-transfer from the electron donor to ZnO 
becomes less favorable, while charge transfer from ZnO 
to the electron acceptor becomes beneficial. As Fig.Qa) 
shows, this results in a turning point at an 0-H coverage 
of approx. 18 %, where the adsorption of hydrogen on a 
zinc atom next to an 0-H bond becomes energetically fa¬ 
vorable compared to the formation of another 0-H bond. 
This argument implies that there is a critical work func¬ 
tion at which the adsorption of H on O and Zn is nearly 
isoenergetic and, thus, there is no distinct energetic ben¬ 
efit of one bond over the other. 

This work function will essentially stabilize itself, as 
the adsorption of H on O lowers the work function and 
facilitates further adsorption on Zn, while conversely, the 
adsorption of H on Zn increases the work function and 
leads to further adsorption of H on O. It would thus ex¬ 
plain why the work-function saturates before the hydro¬ 
gen layer is completed. 

Hence, for coverages around and above 18 %, always a 
complex mixed monolayer with 0-H and Zn-H moieties 
will form. Therefore, it should be emphasized that for 
higher hydrogen coverages, the work-functions calculated 
with the artificial adsorption scheme adopted here do not 
correspond to adsorption structures that will occur in 
experiment. 

Figure [^b) shows the calculated hydrogen-induced 
work function change A4> for a number of different ad¬ 
sorption geometries. To compare the effect of 0-H vs. 
Zn-H bond formation on surfaces with different amounts 
of 0-H bonds already present, A4> is plotted for an 
adsorption pattern where one additional H atom is ad¬ 
sorbed on either another surface O site or a first surface 
Zn site (red diamonds and blue circles, respectively). For 
the (4 X 4)-supercell this means that the work function 
change for a certain 0-H pre-coverage corresponds to an 
overall hydrogen coverage which is 6.25 % higher than 
the given 0-H pre-coverage, and it is always referred to 
the work function of the pristine ZnO(1010) surface. In 
the theoretical case of exclusive 0-H bond formation, we 
compute a work function change of up to A4> = —2.5 eV, 
which corresponds to an overall 0-H coverage of 56.25 %. 
However, in experiment this maximum theoretical A4> 
cannot possibly be achieved, because Zn-H bond forma¬ 
tion sets in well before the 0-H coverage can reach such 
high values, as described above. We expect the stabi¬ 
lization of 4> by the formation of a complex adsorption 
pattern of 0-H and Zn-H bonds for coverages around 
and above 18 % (see red dotted line in Fig. [^. This 
stabilization can be illustrated quantitatively by consid¬ 
ering the effect of a first Zn-H bond on the overall work 
function (blue circles). Remarkably, the difference in A4> 
between H adsorption on a surface Zn vs. a surface O is 
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as large as 1 eV in case of a 18 % 0-H pre-coverage, and 
the resulting = —0.95 eV is comparable to the exper¬ 
imentally determined A^max = —0.65 eV. Because of the 
stabilization of ^ by the formation of a complex adsorp¬ 
tion pattern for hydrogen dosages above 18 %, the cor¬ 
responding calculated A$ is expected to be much closer 
to the experimentally determined value. 

Although these calculations and their underlying ad¬ 
sorption model do not allow us to identify the exact 
surface structure and conformation under experimental 
conditions (this would require combining ah-initio ther¬ 
modynamics with statistical methods such as cluster ex¬ 
pansions or kinetic Monte Carlo modeling, which is be¬ 
yond the scope of the present work), our theoretical 
work unambiguously shows that the 0-H monolayer does 
not complete before the Zn-H monolayer starts forming. 
Rather, the results indicate that both adsorption sites 
compete for hydrogen adsorption, which is expected to 
lead to a complex equilibrium between the 0-H and Zn- 
H (sub)monolayers. This suggests that exposure of a 
pristine ZnO(1010)-surface to hydrogen leads, at first, 
to the formation of 0-H-bonds along with a work func¬ 
tion reduction and CAL formation. Further dosing of H 
increasingly favors the formation of Zn-H-bonds, which 
results in a complex adsorption behavior and a mixed 
adsorption pattern, which is reflected in the change of 
slope of the experimentally determined A^ below 100 L. 
It should be noted that the experimentally determined 
rather moderate shift of the work function leading to 
saturation at A^max = —0.65(5) eV [see FigJ^b)] refers 
to the first ^ini measured after surface preparation. As 
this procedure involves slow cooling of the sample, the 
initial surface H coverage cannot be considered to be ab¬ 
solutely zero. We thus expect ^ini < ^pristine and hence 
A^exp < A^theo, which is in accordance with the exper¬ 
imentally and theoretically determined values. 

As mentioned above, the properties and the function¬ 
ality of interfaces between ZnO and molecular adlayers 
may strongly depend on the local, microscopic electronic 
structure. A complete picture of the H-induced changes 
at the ZnO(1010) surface therefore requires a microscopic 
description of the changes to the electrostatic potential 
[/. As a first step, we calculate the change of U as 

AU{x,y,z) = U^’^^+^{x,y,z)- 

{U'^’^‘^{x,y,z)+ U^{x,y,z)) (3) 

where is the potential of the ZnO surface with 

adsorbed hydrogen, the surface without adsorbed 

hydrogen, and the potential of a hypothetical, free¬ 
standing hydrogen layer in the same geometry as the ad¬ 
sorbed hydrogen atoms. The potential was calculated on 
a regular grid with a grid spacing of 0.1 A in each direc¬ 
tion. The resulting AU is shown in Fig. |^a) along the 
surface normal for two different adsorption geometries: 
First, OH with a single hydrogen in a 4 x 4 supercell 
that is adsorbed on an oxygen atom [corresponding to 
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FIG. 5. (a) Plane-averaged potential change, at the 

ZnO (1010) surface for two different coverage geometries: a 
single 0-H bond [red curve, (b)] and 3 0-H bonds and a 
single Zn-H bond [violet curve, (c)] per 4x4 supercell. The 
location of adsorbed H is highlighted by markers in (b), (c). 


a coverage of 6.25%, (red curve)], and second, OHsZnH 
which contains 3 0-H and one Zn-H bond (violet curve). 
The corresponding adsorption geometries are depicted in 
Fig. [^b) and (c). The latter corresponds to a first Zn-H 
bond on a surface with 18.75% 0-H pre-coverage and it 
is thus consistent with our previous calculations. For the 
sake of simplicity, here we averaged the potential in the 
X and y directions.!^ 

Notably, in both cases the potential changes extend 
less than 20 A into the ZnO bulk, which makes them 
strongly localized along the surface normal compared to 
surface charge accumulation layer de pths i n conventional 
semiconductors such as GaAs or The OHsZnH 

geometry exhibits a stronger confinement and the band 
bending extends only across the first four ZnO layers. 
To our knowledge, this is the first DFT calculation of the 
strongly confined surface downward band bending in such 
a material. It, thus, represents the first quantitative and 
microscopic confirmation of the established qualitative 
view of the strong charge accumulation layer in n-type 
materials [cf. Fig. ma)].!^ 

Apart from the laterally averaged changes to the 
surface potential, we computed the three-dimensional 
changes of U with respect to the pristine ZnO (1010) 
surface for the adsorption geometries displayed in 
Fig. [§b, c). Hence, our DFT calculations also pro¬ 
vide a detailed description of the lateral structure of 
H-induced modifications to the surface potential energy 
landscape. Figures [^b) and (c) show the H-induced po¬ 
tential changes for the OHsZnH adsorption geometry in 
the (x-^)-plane at 2 : = 0 A and in the (x- 2 :)-plane, re- 
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FIG. 6. Hydrogen adsorption induced potential change, At/, 
for OHsZnHi adsorption geometry, (a) Dependence along the 
red dotted line indicated in (b) for different distances from the 
surface. = 0 A corresponds to the line shown in (c). (b) 
AU in false colors along the (x-y)-plane at z = 0 A. (c) AU 
shown for the (x-z)-plane cut along the line indicated in (b). 
Isolines are drawn at 0.5 eV increments of At/. 


spectively. It is shown that 0-H bonds lead to a pro¬ 
nounced reduction of U (blue) which laterally extends 
over roughly 5 nm. This behavior is nearly identical to 
the AU caused by the lone 0-H bond in the OH adsorp¬ 
tion geometry (not shown). The Zn-H-induced increase 
of the potential (red) exhibits a similarly localized char¬ 
acter. As a consequence, At/ is close to zero everywhere 
else. 

Figure a) shows AU along a line that is crossing an 


0-H site and a neighboring pair of 0-H and Zn-H ad¬ 
sorption sites for different distances z from the surface, 
as indicated in the top and side view [Fig. ||:b,c)]. The 
definition of ^ = 0 A is shown in (c). It should be noted 
that z = 0 A corresponds to a position slightly above 
the final ZnO layer and thus roughly agrees with the po¬ 
sition of further possible adsorbates. At ^ = 0 A, AU 
is strongly corrugated, varying between —3.58 eV at the 
single 0-H bond and +1.54 eV at the Zn-H bond. Again, 
the counteracting character of 0-H and Zn-H bonds be¬ 
comes directly visible. As expected, for an excess of 0-H 
bonds, and already shown in Fig. the average At/ is 
negative for this geometry. Moving away from the sur¬ 
face, the corrugation of AU is smoothed, which can also 
be seen in Fig. [^c). This smeared out AU is what, at 
even larger distances (30 A), defines the work function 
change Aff>. 

The microscopic perspective on the surface potential 
modifications obtained from our calculations reveals that, 
in the low coverage regime, changes to the surface elec¬ 
tronic structure are strongly localized both laterally and 
along the z-axis. This corroborates our experimental 
finding that the intensity of the CAL peak is increas¬ 
ing upon H adsorption without a distinct change of its 
spectral shape. During the initial stages of H adsorption, 
the predominant formation of 0-H bonds creates local¬ 
ized and independent metallic sites, supposedly by par¬ 
tially filling the conduction band Zn45 states which are 
shifted below The accumulation of these potential 
wells then leads to the experimentally observed increase 
of the CAL peak [cf. Fig. [^b,c)]. 

The spatial extent of these wells, which are confined 
laterally within ^ 5 nm, indicates at which coverages 
the delocalization of charges within the CAL is to be ex¬ 
pected. As can be seen in Fig. at an 0-H coverage of 
18 %, the potential wells partially start to overlap. An 
increase of the 0-H coverage to 25 % and beyond would 
clearly lead to significant overlap and a generally delo¬ 
calized character of the reduced potential. It should be 
noted, however, that at these coverages, the formation of 
Zn-H bonds has already become favorable (cf. Fig. [^, 
which suggests that the occurrence of a laterally delo¬ 
calized CAL roughly coincides with the beginning of a 
mixed adsorption geometry. 

Until now, the experimental characterization of the 
ZnO (1010) surface focused on observables averaged over 
a certain surface area (CAL intensity, A4>), because of 
the finite spot size of the laser. Now we use a local elec¬ 
tronic state as a local probe of the corrugated poten¬ 
tial: The (partly) hydrogen-terminated ZnO (1010) sur¬ 
face favors the formation of subsurface-bound excitons 
(SX), which are localized at surface potential minima 
(see Ref. El). This excitonic state forms within 200 fs 
after above-band-gap-photoexcitation with the 4.19 eV 
pump photons, as schematically shown in Fig. [^a). It is 
then detected by using a time-delayed 4.63 eV probe pho¬ 
ton, which excites the hole-bound electron above the vac¬ 
uum energy. The spectroscopic signature of the SX lies 
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FIG. 7. (a) SX signatures at a pump-probe delay of 1 ps 

for H dosages from 3 to 44 L, modified from Ref. [O] Intensi¬ 
ties are normalized to the number of electrons excited from 
the valence into the conduction band by the pump pulse (cf. 
Fig.^a)). (b)-(c) Strongly simplified scheme of the ZnO sur¬ 
face showing the correlation of H dosage and the resulting 
density of potential wells with SX density for low coverages. 


close to and below as shown in Fig. It is thus 

isoenergetic with the CAL signature [Fig. [^b)]. How¬ 
ever, the two states can be easily distinguished: While 
the CAL signature is an equilibrium property of the H- 
covered ZnO(1010) surface and can, thus, be measured by 
static photoelectron spectroscopy, the SX require prior 
above-band-gap photoexcitation of the sample and are 
observed in the photoinduced changes of the PE spectra 
(see Ref. m for details). Since, in bulk ZnO, the CBM is 
200 meV above [see Fig.j^a)] and exciton binding en¬ 
ergies in ZnO are on the order of 60 meV, we concluded 
in Ref. [9] that the SX must be connected to regions with 
downward surface band bending, i.e., where the CBM is 
located very close to, or below, Ep. In the present calcu¬ 
lations we show that the reduction of the surface poten¬ 
tial is caused by the formation of 0-H bonds and that the 
resulting downward band bending is strongly localized at 
the 0-H sites. This connection between H-induced sur¬ 
face band bending and the SX makes the latter a sensor 
for the strength and degree of the localization of the sur¬ 
face potential changes: An increase of the H coverage, 
i.e., an enhancement of the number of potential minima 
and, therefore, the number of SX sites should, thus, lead 
to an increase of the SX photoemission intensity. 

The change in SX photoemission intensity with the 
amount of adsorbed hydrogen is shown in Fig. IZla). For 
the depicted, comparably low dosages between 3 and 44 L 
we indeed observe a significant SX signal increase. This 
finding is in agreement with the localized and noninter¬ 
acting character of the potential wells created by 0-H 
bond formation, which we found in our calculations (cf. 
Fig. 1^. In the low coverage regime, these wells accu¬ 
mulate without leading to a spatial overlap of the SX 
species, as schematically depicted in Fig. [^b)-(c). In¬ 
terestingly, here the SX remains largely unperturbed by 
significant changes to the macroscopic properties of the 


sample surface caused by H dosing: The work function 
for the presented H dosages changes by = —250 meV 
and at the same time the CAL intensity increases by a 
factor of ~ 10 (cf. Fig. [^. These experimental findings, 
hence, confirm the strongly localized character of the po¬ 
tential change At/, as predicted by our calculations. 

In addition to the careful study of the SX in the low 
H coverage regime, we also examined the limiting case of 
high H dosages where the distance between H adsorption 
sites is greatly reduced and the surface electron density 
is at its maximum due to formation of the CAL. Both 
effects are expected to reduce the SX intensity because 
of (i) spatial overlap and resulting mutual screening of 
the electron-hole pairs and (ii) screening by the increased 
amount of electrons (cf. Refs. [9] and [18]). Indeed, by in¬ 
creasing the H 2 dosage to 150 L, where the CAL intensity 
is at its maximum (see Fig.[^, the SX signal is quenched, 
i.e., we observe barely any pump-induced signal below Ep 
(not shown). 

This observation suggests that in the hydrogen dosage 
regime between 44 L and 150 L, there is a transition of 
AU from the strongly localized character to a laterally 
smeared out CAL which results in the formation of a 
two-dimensional electron gas at the surface as suggested 
by Ozawa et al^ for a comparable coverage. In con¬ 
junction with our theoretical results which suggest an 
enhanced overlap of potential wells for 0-H coverages 
above ~ 25 %, we conclude that the delocalization of 
the surface charge is reached between 44 and 150 L. Re¬ 
markably, this roughly coincides with the coverage regime 
below 100 L at which isoenergetic Zn-H bond formation 
sets in (at 18 % 0-H coverage), as deduced from the work 
function change further above. 

For even higher hydrogen coverages we expect a signif¬ 
icant difference between our theoretical model of adsorp¬ 
tion and the situation in experiment. This is, because 
high hydrogen coverages were shown to lead to drastic 
changes in the surface geometric structure of various ZnO 
surfacethat cannot be accounted for in our model. 
Diffusion of hydrogen into the ZnO bulk may increase the 
number of defects and lead to substantial changes in the 
(near-)surface electronic structure.^^In particular, it was 
found for the Zn-terminated ZnO surface, that hydrogen 
at high dosages forms bonds with initially sub-surface O 
atoms, indicating a complex reconstruction behavior at 
this surface that is suggested as the main stabilization 
mechanism for this polar surface.^^^^ 

Nevertheless, our experimental observation of a re¬ 
duction of the CAL intensity for dosages above 150 L, 
as shown in Fig. is consistent with previous 
experiments.^ Furthermore, our calculations for the 
completely H-covered surface predict, in agreement with 
a previous study,II2l a re-opening of the band gap at the 
surface, which leads to a restoration of the semiconduct¬ 
ing state (see above). Experimentally, however, we never 
observe a complete quenching of the CAL, even for H 2 
exposures up to 500 L. This strongly suggests that, at the 
hydrogen pressures used in our experiments, the hydro- 
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gen coverage of the ZnO(lOlO) surface never completely 
saturates. 


IV. SUMMARY AND CONCLUSION 

A detailed microscopic view of hydrogen adsorbed on 
the ZnO(lOlO) surface has been developed using surface 
sensitive photoemission experiments and DFT calcula¬ 
tions. In the experiment, we observe the emergence of a 
photoelectron signature at the Fermi level which, at first, 
is enhanced by increasing the hydrogen dosage. These 
occupied electronic states are attributed to the gradual 
formation of a charge accumulation layer at the ZnO sur¬ 
face. For dosages exceeding 150 L, the CAL intensity is 
reduced. Simultaneously, we observe a reduction of the 
work function from 4.50(5) eV to about 3.8 eV. We show 
that the observed work function change and CAL inten¬ 
sity variation can be explained by a competition between 
0-H and Zn-H-bonds. Zn-H bonds become energetically 
favored when « 18% of the ZnO(1010) surface is covered 
with 0-H bonds, eventually leading to a complex mix of 
adsorption sites. 

The corresponding changes to the surface potential 
stem from the amphoteric character of hydrogen at the 
O and Zn surface sites. At the former it acts as elec¬ 
tron donor, leading to a reduction of the electrostatic 
surface potential, whereas, in the Zn-H bond, hydrogen 
accepts an electron, thereby increasing the surface po¬ 
tential. The prevalence of 0-H bonds in the low cov¬ 
erage regime causes the crossing of CBM and Ey and, 
thus, the surface metallicity. We calculate that, along 
the surface normal, the electrostatic potential change AU 
reaches only a few nm into the ZnO bulk, which is a di¬ 
rect confirmation of the conventional macroscopic view of 


the charge accumulation layer at the ZnO surface. Lat¬ 
erally, the potential changes in the low coverage regime 
remain mainly confined to the H adsorption site. This 
strong localization of the H-induced AU is confirmed ex¬ 
perimentally by using the signature of the SX as a local 
probe. 

Our study shows that the interaction of hydrogen with 
the ZnO (1010) leads to a complex electrostatic potential 
landscape that has to be treated on a microscopic level. 
The strongly localized character of H-induced changes to 
the surface potential for low H coverages is relevant re¬ 
garding interfacial energy level alignment with functional 
molecules, especially considering that H is a ubiquitous 
contaminant. Thus, when dealing with hybrid interfaces 
of ZnO with functional molecules, not only macroscopic 
surface properties have to be considered, but also the 
electronic structure at the binding site of the adsorbate. 
Knowledge of this effect is therefore crucial for ZnO- 
based device applications and could also be intentionally 
used to tailor the properties of the ZnO surface. 
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